Astro-physical Bulletin, 2008, 63, 389-396 



© Special Astrophysical Observatory of the Russian AS, 2009 



Catalog of Radio Galaxies with z > 0.3. I: Construction 
of the Sample 



o 
o 

> 

o 



O 

U 

Or 

6 



On 
O 



13 



M.L. Khabibullina a , O.V. Verkhodanov a 

Special Astrophysical Observatory of the Russian AS, Nizhnij Arkhyz 369167, Russia 
Received October 2, 2008; accepted October 24, 2008. 



The procedure of the construction of a sample of distant (z > 0.3) radio galaxies using NED, 
SDSS, and CATS databases for further application in statistical tests is described. The sample 
is assumed to be cleaned from objects with quasar properties. Primary statistical analysis of the 
list is performed and the regression dependence of the spectral index on redshift is found. 
PACS: 98.52.Eh, 98.54.-h, 98.54.Cr, 98.62.Ve, 98.70.Dk, 98.80.Es 



1. INTRODUCTION 

Radio galaxies, which are among the most powerful 
cosmic objects, allow us to study the evolution of mat- 
ter and the dynamics of the expansion of the Universe 
at different cosmological epochs. 

Investigations of these radio sources are linked to 
several cosmological tests, which permit estimating 
the parameters and evolutionary characteristics of the 
Universe. Radio galaxies are identified with giant el- 
liptical galaxies with absolute magnitudes M ~ —26 
and have black holes with masses on the of 10 9 Mq 
as their central engines. These facts allow the radio 
galaxies to be used as tools for the study of the pa- 
rameters of the distribution of visible and dark mat- 
ter, dynamics of the Universe, and the history of the 
formation of its structure. Let us point out various 
published methods developed to find cosmological pa- 
rameters from the observed characteristics of radio 
galaxies [TJ [2j [3] , such as: 

• "magnitude-redshift" Hubble diagrams 
(«K-*») 0J; 

• the "size-z" diagram 0[6l[7]; 

• the "age-z" diagram [5]; 

• the "source number-flux' (log N — log S) dis- 
tribution [91 [TO]: 

• search for clustering and determination of the 
clustering parameter [TT] ; 

• search for gravitational lenses [12] ; 

• investigation of the properties of the radio 
haloes of clusters of galaxies [13] ; 

• cosmic background temperature variations 
caused by the interaction of the large-scale structure 
(clusters of galaxies) at different redshifts with cosmic 
microwave background and observed on small angular 



scales (the Sunyaev-Zel'dovich effect due to inverse 
Compton scattering) [14] and on large scales in vari- 
able gravitational potential — the Sachs-Wolfe effect 
[15] . Radio galaxies can be used as a tool for search- 
ing for clusters of galaxies [TO] : 

• and a number of other methods (see, e.g., re- 
views Huang). 

A striking feature of radio galaxies is that we ob- 
serve them virtually from their formation epoch, i.e., 
their radio sources are so powerful that modern ra- 
dio surveys have catalogued almost all objects of this 
type [3]. Such objects can serve as good probes for 
the study of the formation of clusters of galaxies. 
Thus Venemans et al. [TB] show that 75% of the ra- 
dio galaxies with z > 2 are associated with proto- 
clusters. Based on this fact, the above authors esti- 
mate that there arc roughly about 3x 10~ 8 of forming 
clusters per comoving Mpc 3 in the redshift interval 
2 < z < 5.2 with an active radio source. They also 
estimate the density of protoclusters within the same 
redshift interval to be of the order of 6xl0 -6 Mpc -3 . 

Identification of radio galaxies with giant elliptical 
galaxies (gE) born as a result of mergers in the early 
epoch allows them to be used not only to probe the 
formation of the large-scale structure and test mod- 
els of star formation. The ( "K-z" ) diagram confirms 
the age of the old stellar population in gE-type galax- 
ies [4] and this fact can also be used to compute the 
ages of these galaxies [HI HU [20] and perform quick 
photometric redshift estimates for these objects out 
to z ~ 4 [2"T] . 

Another important problem associated with the 
study of radio galaxies is the origin of supermas- 
sive black holes. The estimates of energy release 
in the central engines of radio galaxies at redshifts 
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4 < z < 5, when the Universe was slightly older than 
one billion years (in terms of the ACDM model) , yield 
black-hole masses on the order of 10 9 Mq, whereas 
it takes much more than one billion years for such 
black holes to form [22]. This paradox can be re- 
solved in terms of the hierarchical formation model, 
but the actual number of such objects (black holes 
with masses on the order of 10 9 Mq in the redshift 
interval 4 < z < 5) in the Universe is not known with 
certainty. 

Distant radio galaxies are selected by a criterion 
based on the spectral index of their continuum. To 
this end, lists of sources with ultrasteep spectra are 
compiled, where the parameters needed are the low- 
frequency data points on the continuum radio spec- 
tra. Among such objects a high percentage of distant 
radio galaxies is found [H [H [11 HI US] , including 
the most distant objects with redshifts z > 4.5: 
z = 5.199 [22] and with z = 4.514 [28] . 

Nowadays databases and identification tech- 
niques allow automating search and identifica- 
tion of distant radio galaxies. A number of au- 
thors [29] [301 El O [331 (311 [35] applied such ap- 
proaches to study radio sources using CATS0 [51 [37] 
and NED El (NASA/IPAC Extragalactic Database) 
databases of astronomical catalogs and Sky View vir- 
tual telescope d, which allow performing morphologi- 
cal analysis of objects. 

The catalog of distant radio galaxies can be used 
not only to perform just purely cosmological stud- 
ies, but also carry out statistical analyses of iden- 
tification lists and of the corresponding populations 
at different wavelengths [381 El [H [41] ; to search for 
and analyze the properties of subsamples of radio 
galaxies (42J, [43l HH [45l [35], and simulate radio- 
astronomical surveys at RATAN-600 [H S3 [48] . Or- 
ganizing the databases and identifying the available 
heterogeneous measurements with particular objects 
are also tasks of great interest, and, as practical ex- 
perience shows, such compiled catalogs are popular 
among the astronomical community (37j [H [50] . 

Moreover, radio galaxies, as extended ob- 
jects, may distort microwave background and bias 
the estimates of the angular power spectrum 
of CMB [53(33]. 

In this paper we report a catalog of radio galax- 
ies with spectroscopic redshifts z > 0.3 based on the 
data from the NED, SDSS3[52 , CATS, and other 
databases and on an analysis of the corresponding 
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bibliography, which allowed us to discard objects ex- 
hibiting quasar properties. There is no yet a consen- 
sus among the research community as to which radio 
galaxies should be considered "distant", with differ- 
ent authors using this term to denote objects with 
redshifts above 0.3, 0.5, or 1.0. Based on observa- 
tional data, one can set the boundary for distant ra- 
dio galaxies at z ~ 0.6 — 0.7 (corresponding the age 
of the Universe of about 6.5 Gyr), when, according 
to ACDM cosmology, the domination of dark energy 
(DE) began. The lower redshift boundary can then 
be set at half this z value. Such a choice permits the 
test to include the population of galaxies after the 
beginning of the DE domination. 

This is the first publication of a series of three 
papers dedicated to the construction of the sam- 
ple and statistical analysis of the catalog of radio 
galaxies, which contains objects with known spectro- 
scopic z and available radio data and optical photom- 
etry. In the future we plan to perform cosmological 
tests based on the physical parameters of the objects 
of our catalog. 

2. THE CATALOG 

2.1. Selection of Objects 

To build our primary list, we use the NED database, 
where we select objects with the following parame- 
ters: redshift z > 0.3 and morphological type "radio 
galaxy" . The initial list contained a total of 3364 
objects. This galaxy sample is polluted by objects 
with incomplete data or objects with different 
properties. Our next stage therefore consisted in 
cleaning the initial sample from wrong sources. 
To this end, we select the following objects to be 
removed from the initial list: (1) objects with pho- 
tometrically determined redshifts; (2) objects with 
quasar properties. We performed an extensive scan 
of the literature including optical, infrared, and radio 
surveys represented both in the NED and CATS 

[51[51[210230iDG^ 

El M El EH El El El El EJ EtJ EJ EJ [HHUm 

us eh mi m m mi na ed isaEiiMiEsiEiEa 

i^iMimiiTmiiTmiirmi^ 

HH HIU HH Q331 QT4]. The complete catalog is avail- 
able from the CATS database!!, and CDS ftfl CDS 
databasfl Our sample also includes radio galaxies 
investigated within the framework of the 
"Big Trio" [2TH 1271 11T51 im [IT71 im 
program, i.e., objects of the RC catalog [121[ 1122] 



5 |f tp : //cats . sao . r u/pub/CATS/RGLISf] 

6 ftp : //cdsarc . u-strasbg . f r/pub/cats/ J/other/AstBu/64 . 123 

7 http : //cdsarc . u-strasbg . f r/cgi-bin/VizieR?-source=J/other/AstBu/64 . 123 



CATALOG OF RADIO GALAXIES WITH z > 0.3. I: CONSTRUCTION OF THE SAMPLE 



391 



Figure 1: Sky positions of the selected radio sources in Galactic coordinates. The white circles and gray crosses 
indicate the SDSS objects and all other sources, respectively. 



with published rcdshifts and measured photometric 
magnitudes. 

As a result, we compiled a catalog containing 2442 
objects for further studies. 

Figure 1 shows the location of selected radio 
galaxies on the sky. 

2.2. Description of the Catalog 

The catalog of radio galaxies consists of three parts 
with each dedicated a separate paper. Here we give 
the list of the objects with the names of the source 
catalogs (and, in the case of the 3C and 4C sur- 
veys, with the source names); equatorial coordinates 
(J2000.0); spectral indices for the 325, 1400, and 
4850 MHz frequencies, and spectroscopic redshifts. 
No spectral indices are given for galaxies with the 
fluxes measured only at a single frequency (e.g., SDSS 
objects with the data from the FIRST catalog [84]). 
We compute the spectral indices based on the results 
of cross identification in the CATS database with 
a 200" x 200" identification box. To remove spuri- 
ous field radio sources within the given box, we ap- 
ply the technique of data analysis similar to that 
described by Verkhodanov et al. [221 [51] . 

The central idea of the method consists in per- 
forming joint analysis of the data in the coordinate 
and spectral spaces in order to identify likely iden- 
tifications of particular radio sources at different ra- 
dio frequencies. To this end, we use the 'spg 1 pro- 
gram [123J of the system of the reduction of contin- 
uum data at RATAN-600 radio telescope. To repre- 
sent the spectra (S(v)) for the subsequent computa- 
tion of spectral indices we parameterize them by the 



following formula lg S{v) = A + Bx + Cf(x), where 
S is the flux in Jy; x, the logarithm of frequency v 
in MHz, and /(x) is one of the functions exp(— x), 
exp(x), or x 2 . 

We collected the optical (DSS 2) and radio im- 
ages (mostly from the NVSS [96]) of the galaxies 
and their spectra into an atlas, which is available at 
http : / / sed . sao . ru/ ~rita/Atlas_RG . html. 

2.3. Statistical Analysis of the Sample 

Figure 1 shows the locations of the selected galax- 
ies in the sky. The list as a whole is by no means 
complete or homogeneous, since the sample contains 
objects brought from different catalogs and located 
in mutually unconnected sky areas. The most com- 
plete and homogeneous subsample of our list contains 
SDSS objects, which are indicated by white circles in 
Fig. 1. Among the SDSS radio galaxies (with the data 
from the NVSS and FIRST radio surveys) there are 
many objects with small redshifts (z < 0.5) and small 
fluxes (S < 15mJy), making the subsample in ques- 
tion stand out qualitatively against other subsamplcs. 

We nevertheless view the list reported here as a 
tentatively pure catalog of radio galaxies suitable for 
testing cosmological relations. We call this list only 
tentatively pure, because we cannot say with cer- 
tainty that selected objects would not be found to 
exhibit properties typical of another class. The list 
remains open with both some objects to be possibly 
removed or new objects added. 

The Table and Figs. 2-6 give the primary statis- 
tical parameters for the resulting catalog. 

Figure 2 shows the histogram of spectral indices at 
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Table 1: Table. Parameters of the sample of the first 
part of the catalog of distant radio galaxies (z > 0.3): 
the median reshift, minimal, median, and maximum 
1400 MHz flux in mjy, the minimum, median, and 
maximum spectral index at 1^00 MHz 
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Figure 2: Histogram of spectral indices at 14-00 MHz. 

1400 MHz. We fitted the spectra of the sample mostly 
by linear functions, resulting on almost wavelength 
independent distribution. The median of the distri- 
bution of the spectral indices for the entire sample 
is equal to -0.63. Figure 3 shows the distribution of 
fluxes at the same frequencies H, and Fig. 4 shows the 
distribution of redshifts. 

Figures 5 and 6 show the "flux-redshift" and 
"spectral index-redshift" diagrams, respectively. 

The "flux-redshift" relation shows a well-defined 
upper flux envelope (at 1400 MHz, NVSS data), 
which gives the maximum luminosities of observed 
radio galaxies at different z, and apparently is due 
to the dynamics of cosmological expansion. We fit 
linear regression to the maximum values of the dis- 
tribution in Az = 0.5 wide bins. Our choice of the 
Az = 0.5 bin width was determined by the following 
reasons: it must be (1) small enough (Az < 1.0) to 
take the differential properties of the population into 
account and (2) big enough (Az > 0.1) for the bins 
to contain sufficient number of objects. We adopted 



8 We will report the fluxes in the third paper of this series. 
The plots here are based on the data adopted from the NVSS 
survey. 
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Figure 3: Histogram of the fluxes of radio galaxies 
according to NVSS data. 
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Figure 4: Histogram of redshifts. 

Az = 0.5, which is, on the one hand, halfway be- 
tween the two extremes, and, on the other hand, en- 
sures that bins at z > 3 contain a comparatively 
large number of objects. The regression relation is 
given by the formula log Sunn = p + rz, where S1400 
is the 1400 MHz flux; p — 4.16 is a constant intercept, 
and r = —0.45 is the slope of the line. The objects of 
this sample are to be divided by morphological type 
and luminosity for further analysis. Here we point out 
the especially powerful radio galaxies with fluxes lie 
by more than half a magnitude above the regression 
boundary: 3C 295 (z=0.464) [102], PKS 0742+10 
(z=2.624) PHI, 8C 1435+635 (2=4.261) [315]. We 
will analyze the properties of these objects in the 
third paper dedicated to the catalog. We also do 
not plot the data point for the radio galaxy VLA 
J123642+621331 (z=4.424) 126], which is included 
into the catalog, but is absent in the NVSS list (be- 
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Figure 5: TTie "observed flux-redshift" diagram for 1400 MHz (the NVSS catalog). The triangles and circles 
indicate the SDSS data and objects of all other catalogs, respectively. The regression fit is drawn for the 
maximum fluxes inside the bins with a redshift, step of Az = 0.5. 



cause of its weak flux (0.5 mJy)). 

The "spectral type-redshift" diagram (Fig. 6) 
shows a well-defined trend — spectral index de- 
creases with increasing z and fits linear regression 
relation a = a + bz, where a = —0.73 ± 0.02 and 
b = —0.15 ± 0.01 are the intercept and slope, respec- 
tively. Some objects exhibit peaks near 1 GHz and at 
lower frequencies, which result in positive spectral in- 
dices at low frequencies. In cases of negative a, with 
steep source spectra, we fitted the spectrum by a lin- 
ear relation. The spectral indices in the diagram is 
computed at 1400 MHz. Selection of objects by the 
steepness of their spectra has been used extensively 
by a number of authors [24], [26l [27l [25] , but we are 
the first to derive an analytical regression formula for 
such a large sample of radio galaxies. We compute 
the regression relation for the median of spectral in- 
dices in z bins with a bin width of Az — 0.5, which 
we adopt based on the same considerations as in case 
of "flux-redshift" diagram (see above), thereby mak- 
ing the inferred regression parameters stable against 
the eventual future enlargement of the sample. The 
slope can, in principle, be explained by two selection 
effects, although physical causes cannot be ruled out: 

(1) We adopt the lists of distant radio galaxies 



from the catalogs compiled by selecting objects with 
steep spectra. The only way to avoid this effect is to 
use complete samples of radio galaxies with measured 
rcdshifts, which are few. 

(2) When we select the most powerful radio 
sources at large z the radio flux is contributed by 
hot spots, which have steep spectra. And the farther 
is the radio source, the more likely for it to be more 
powerful compared to the surrounding sources and to 
have a steep radio spectrum. We plan to verify this 
fact in our future studies by analyzing subsamples 
of radio galaxies of different luminosities. However, 
despite the scatter of spectral indices, the analytical 
form of the regression can be used for preliminary 
selection of radio galaxies and obtaining preliminary 
distance estimates, the regression yields a median es- 
timate for a large number of objects. 

3. CONCLUSIONS 

We report the first part of the catalog of radio galax- 
ies with z > 0.3 with complete information on the 
coordinates of the objects, their spectral indices at 
radio frequencies, and redshifts. We constructed our 
catalog based on the lists of objects from the NED 
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Figure 6: The "spectral index-redshift" diagram (the spectral indices are computed at 14-00 MHz). The triangles 
and circles indicate the SDSS data and objects of all other catalogs, respectively. The regression fit is drawn 
for the median spectral indices inside the bins with a redshift bin size of Az = 0.5. 



and CATS databases. Our sample includes objects 
of the SDSS survey and of the "Big Trio" program. 
We perform preliminary statistical analysis and com- 
pute the distributions of spectral indices and fluxes 
and plot diagrams of parameters as functions of z. 
The catalog forms a basis for further study of the ob- 
jects and of the properties of distant radio galaxies 
population. For our objects we construct the spectra 
and the "flux-redshift" and "spectral index-redshift" 
diagrams. On the former diagram we determine the 
upper boundary of the dependence of flux on z. 
which can be described by the regression relation 
log S = 4.16 — 0.45z. The latter diagram yields the 
regression relation a(z) = —0.73 — 0.15z. 
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